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Abstract 

^\ 

■ We present, solve and numerically simulate a simple model that describes the con- 
, sequences of increased longevity on fertility rates, population growth and the distri- 
ct \ bution of wealth in developed societies. We look at the consequences of the repeated 

use of life extension techniques and show that they represent a novel commodity 
whose introduction will profoundly influence key aspects of economy and society in 
I \ general. In particular, we uncover two phases within our simplified model, labeled 

^ ' as 'mortal' and 'immortal'. Within the life extension scenario it is possible to have 

. sustainable economic growth in a population of stable size, as a result of dynamical 

\ equilibrium between the two phases. 
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1 Introduction 



Constant lifespan is often taken as one of the assumptions in creating and an- 
alyzing econophysical and sociophysical models. Today's bio-medical research 
is, however, uncovering the reasons why organisms live as long as they do. 
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Through modem genetic engineering, the apphcations of this research are 
converging on the point of finding practical ways to extend hfe substantially 
(and possibly repeatedly) beyond current life expectancy. 

The uncovering of the so called "secret of life" was one of the crowning achieve- 
ments of the second half of the past century. The discovery of the structure 
of DNA by Crick and Watson [1], and the later successful translation of the 
DNA code into the language of proteins, fueled the continuing revolution in 
molecular biology and bio-technology. This revolution is now making it pos- 
sible to tackle rationally the complementary question of why we die. For the 
first time we have the option of looking at death from a fact-based [2] sci- 
entific perspective. The picture that is emerging is quite unexpected [3]. We 
are getting farther from the concept of "natural death" as an immutable and 
inevitable end of life 

Physical models played a crucial role in the discoveries that marked the birth 
of molecular biology [7]. Similarly, the well-developed machinery for under- 
standing the behavior of complex systems is today being positioned to help 
understand the mechanisms of death and the extension of life [51I^1U|11] . An- 
other avenue of research is to try to understand the social and economic impli- 
cations of the prolongation of human lifespan. In principle, there are two quite 
distinct paths that one may choose to take [12]: economics or econophysics. 
In this letter we approach these problems from an econophysical framework. 
Using this language we show that it is possible to model and predict some of 
the far-reaching social and economic consequences of the successful extension 
of human lifespan that have, up to now, been disregarded both by economists 
p^Sfl^ and physicists p^l6|17] . 

The fact that we all must die has been one of the central points shaping all 
human societies. Substantial modification or even the removal of this mor- 
tality paradigm will necessarily bring about great change in how societies 
function. It is important to try to anticipate these changes. Successful model- 
ing of these phenomena is not only of practical, but of heuristic value. Many 
important discoveries, particularly in physics, have followed from analyzing 
the consequences of modifying key paradigms. The introduction of a new and 
extremely sought-after commodity, allowing for the extension of life, would 
bring about a great change in economy and society in general. 

In this paper we present and solve a simple model dealing with the conse- 
quences of just such a novel commodity. We study the implications of possible 
long-term extensions of life on society and its economy. We model the dy- 
namics of social and economic indicators of a society and investigate how the 
introduction of life extension will influence fertility rates, population growth 
and the distribution of wealth. For this purpose we propose and analytically 
solve a simple model. The presented model, when life extension is absent, is 
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related to earlier investigations [T8p9] . the main new feature being the in- 

troduction of overall economic and population growth. The presented results 
include conclusions that the population explosion is not a necessary conse- 
quence of introduction of life extension commodity, and that it is even possible 
to have sustainable economic growth in a population of stable size. 



2 Basic model 

Our model tracks through the generations the number and individual wealth of 
all of the descendants of a specific individual at generation t = having wealth 
m. We begin by first treating the simpler case of no life extension. Within our 
model, the life span of each individual consists of three phases: formative years 
(parents invest in the individual); adult years (individual inherits some initial 
money, marries, the pair earns some final amount of money, has children); 
old age (individual lives off his pension and ultimately dies). As a result of 
these assumptions, dependency ratios are fixed within our basic model, and 
are constant throughout the population. We track the adult phase of each 
individual which starts at t and lasts until t + 1. We assume that individuals 
inheriting m money choose spouses having the same amount of money, i.e. that 
the pair starts off with 2m initial capital. While this may appear to be a natural 
assumption it is not obvious if it holds empirically. For example, Dragulescu 
and Yakovenko [20j have studied the related phenomena of the earnings of 
spouses and have shown the earnings to be essentially uncorrelated. It would be 
interesting to investigate the correlation of inherited wealth of spouses. Within 
the model presented here we stay with the above simplifying assumption. We 
further assume that society is numerous enough so that everyone can find a 
mate. During their working life the pair increase their wealth by a factor 7, a 
fixed constant for the whole society satisfying 7 > 1. This money is spent on 
their children and the pair's pensions: 



In the above equation k is the number of children, C the investment in each 
child, m' the inheritance of each child. The pension is assumed to be propor- 
tional to initial wealth, i.e. P{m) = am. The number of children k is assumed 
to take on the maximal possible value consistent with the rule m' > m. This 
is a crucial assumption strongly affecting the model's predictions. It implies 
that parents have children only if they can assure them equal or better finan- 
cial start-up compared to what they had. The number of children thus follows 
from a simple economic criterion. As a consequence, the model leads to a 
positive relation between fertility and wealth. For this reason, it is obviously 
not applicable to poor societies. In those societies the choice of the number 
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Fig. 1. (a) Numerical simulations of the time evolution of wealth distributions for 
the case of no life extension. The plot corresponds to 7 = 2, a = 0, C = 2, and an 
initial Gaussian wealth distribution, (b) Comparison of the time dependance of the 
fertility rate for France ^33j with the predictions of the model without life extension 
when 7 — a decreases slowly (adiabatically) bellow the critical half-integer value. 

of children is more strongly related to survival (procreation) and less to ex- 
pectations of their future wealth. Our ultimate goal is to analyze the effects 
of life extension on society. For this reason we focus on developed societies in 
which economic choices play a dominant role, societies having the necessary 
financial means to purchase life extension. Although many other things influ- 
ence fertility (e.g. religion, level of education, system of beliefs), the criterion 
chosen is that of a simplified model seeking to capture the dominant aspect 
of the relation between fertility and wealth in developed societies. 

We will keep track of m{t) and n{t) (money inherited by the descendants at 
generation t and the number of those descendants). Note that the total number 
of people in the society at time t is simply N{t) oc n{t)/2*. 

It is easy to see that for non-trivial dynamics we need to further have 7 > 
a + 1/2 + C /2m, since smaller values of 7 lead to = for all values of m. 
By introducing the critical value m^, = C/(27 — 2a — 1), as well as auxiliary 



quantities Ki = ^"^(j ") and Mi = ^'"'^^ — C, we can write the solution 



C+m 



of the above dynamics as 



2m{-j—a) 



k = 6{m — m^,) Ki , 
m' = 6{m — m^,) Mi . 

Square brackets denote the integer part of an expression. The step function 
6{x) used here equals unity for x > 0, and vanishes for x < 0. 

The above equations make it possible to investigate the dynamical evolution 
of wealth distributions from given initial conditions. Wealth distributions have 
been extensively studied in the literature. The field began with the power law 
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distribution of Pareto [2T]. Recent investigations show that, while Pareto's law 
gives a good fit at higher incomes, it does not agree well with observed data at 
middle and low incomes P^23f24ll25f26] which best fit to lognormal or Gibbs 
distributions. The two regimes follow ^7112^ from the fact that in low and 
middle income ranges the accumulation of wealth is additive, while in the high- 
income range wealth grows multiplicatively. Simplified models [28ll29|30j have 
linked Gaussian wealth distributions with egalitarian societies. Our model 
without life extension agrees well with this phenomenology. The dynamics of 
the model is such that it preserves Gaussian shaped wealth distributions, as 
shown in Fig. [T^. 

In contradistinction to this, Sala-i-Martin has shown [3T][32] that highly segre- 
gated societies exhibit bi-modal wealth distributions. In the next section, we 
will show that the introduction of life extension in our model can lead to so- 
cietal segregation, which then results in the appearance of just such bi-modal 
distributions. 

Within the framework of our model it is also possible to analyze fertility rates 
of a given society. The fertility rate / is the average of k{m) over the whole 
population. If the wealth distribution is such that most of the population have 
wealth C then the above solution gives / < [2(7 — a)]. Due to the integer 
part operation, the fertility can depend strongly on small changes in economic 
growth 7 or of social expenditures a. If 7 and a change slowly with time, 
then 7 and a can be neglected it the equations of motion and we uncover the 
same relation between fertility, economic growth and social expenditure. As a 
result, even the smallest decrease of 7 — a bellow half integer values leads to a 
decrease of fertility by one unit. This is illustrated in Fig. [T]d. The data points 
correspond to measured fertility rates in France. Similar abrupt decreases of 
fertility have been observed for many other developed countries [33j- Thus, an 
increase in social expenditures greater than the increase in economic growth 
results in a step-down in fertility rate. The sharp decline of fertility rates in 
developed countries has been most often accredited to increased participation 
of women in the work market. The study of the relation between fertility and 
wealth within the presented simplified model may offer new insight into this 
important phenomenology. 

It is important to note that the simplifying assumptions made in this paper 
make the presented toy model (and its generalization to include life extension) 
analytically solvable. Future models will need to be made more realistic. To do 
this it will be necessarily to relax some of the assumptions of the current model. 
In particular they will have to treat the effects of overlapping generations. 
Unlike in our mean-field model, a more realistic model will have to have agents 
with different growth factors, different life spans, richer or poorer spouses, etc. 
We intend to study the effects of the relaxation of these assumptions in a future 
publication. These more realistic models will have more phenomenological 
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input parameters and will necessitate a purely numerical treatment. We hope 
that the present analytically solvable model will serve as a useful zeroth-order 
approximation to these future models. 



3 Model with recursive life extension 



We now generalize the model to include life extension. The extension of life 
for one individual and one time step (equal to the natural length of the adult 
period) costs E. We assume that any individual having enough money to pay 
for this life extension will do so, no matter what. Therefore, for 7m > E we 
now have: 

27m = 2E + kC + km + 2m' . (2) 



On the right hand side the first term pays for life extension for the pair, the 
second and third terms are the investment and inheritance of each child, while 
the last term represents the "inheritance" of the original pair with which they 
begin their new life cycle. Note that we have assumed that both parents inherit 
the same money as do each of their children. After life extension individuals 
are assumed to be working able, i.e. there is no pension term in this case. 
Note that for 7m < E we have the same dynamics as before, i.e. as given 
in Eq. ([1]), with P{m) = am. Note that the introduction of life extension 
decreases dependency ratios - in the extended life periods the population is 
of good health and is assumed to be in the economically productive phase. 
Recursive application of life extension drives the dependency ratio to zero. 

We introduce mi = E/'y and m2 = E/{'y — 1). The life extension phase is for 
m > mi. The solution of the model depends on the relation between critical 
values and mi. In terms of K2 = [ '"^^~^^~'" ] and M2 = ^""^'j^f-/'^' , the 
solution for < mi is given by: 



k = {6{m — m^) — 9{m — mi)} Ki + 9{m — m2) K2 , 
m' = {6{m — m*) — 6{m — mi)} Mi + 6{m — 7/12) M2 . 

Similarly, the solution for m^ > mi equals: 



k = 6{m — max{m2, m^,}) K2 , 
m' = 6{m — mi) M2 . 

Note that for all values of m the function 6{m — mi) measures if life was 
extended in the current generation. 
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From the above solutions it follows that, in the life extension phase, A; 7^ is 
possible only for the case 7 > 3/2, and m > (C + 2E)/{2'y — 3). We also see 
that m' > m everywhere except for m G [mi, 7712). It follows that individuals 
with m > 1712 are immortal - they extend their lives and later they (and their 
children if they have any) have more money then in the previous life cycle. 
Note that the potential segregation between mortals and immortals can lead 
to serious political tensions and instabilities in a society. We next tackle the 
issue of segregation. 

Fig. [2] illustrates the obtained solutions for the time dependence of the number 
of children k and wealth m for a society with 7 = 3/2, a = 1/2, C = 1. Fig. [2^ 
and [2)d correspond to the case of no life extension. The part of the population 
with m < 1 has no children and dies off, those with m > 1 have one child. 
That child is financially better off than its parents, i.e. m' is above the dashed 
m' = m line. After t generations we have m{t) ~ 2*, while the population 
decreases as N{t) ~ 2~*. Fig. ^ and correspond to the same society as 
before but with life extension costing E = ?>. The population now consists of 
two groups that never mix - mortals and immortals. Immortals have m > 6. 
The number of mortals roughly decreases as 2~* while their individual wealth 
oscillates in the interval [1,6). In fact, numerical simulations show that the 
majority of this population oscillates in a narrow interval of wealth around 
m = 2. For k = Q (immortals without children) we have m' = 3m/2 — 3, 
so that their wealth grows asymptotically as (3/2)*. For > the condition 
m' > m is not met, i.e. for the case considered immortals can have no children. 
From Eq.([2]) it follows that immortals can procreate only if economic growth 
is such that 7 > 3/2. 

Even this single example shows how the introduction of life extension severely 
affects both population growth and wealth distribution of the whole society, 
not just of the newly created class of immortals. The key effect of life extension 
on mortals is that for some of them life span is increased - some do cross 
into the life extension phase m > 2, however, their wealth then decreases 
(solid curve bellow the dashed m' = m line) and ultimately makes further life 
extension impossible. 

Fig. [3] illustrates the time dependence of wealth for two different sets of pa- 
rameters. The interval [mi, 7712) represents a barrier through which a mortal 
must "tunnel" in a single generation in order to become immortal. The only 
way that the descendants of mortals can become immortal is if m = mi — 0+ 
leads to m' > m2. To get this we must have: 



m* < mi 



2mi(7 ~ Q;) 
C + m2 



> 



2mi(7 



a 



C + mi 
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Fig. 2. Solutions of the model for 7 = 3/2, a = 1/2, C = 1. Plots (a) and (b) 
correspond to no life extension, plots (c) and (d) to life extension with E = 2>. The 
dashed lines in (b) and (d) correspond to ml = m. Note that the introduction of 
life extension has decreased fertility. 



(a) 



(b) 




Fig. 3. Time evolution of individual wealth: (a) An example of a society where 
tunneling is not possible and mortals and immortals remain segregated (7 = 3/2, 
1/2, C = 1, E = 3); (b) An example with tunneling (7 = 2, a = 0, C = 2, 



a 



E = 3). Shaded regions denote corresponding [mi, 7712) intervals. Horizontal lines 
beneath them denote the critical value m = m^. 

Equivalently, these inequalities may be written as: 



2§(7-«) 
' ^ C7-I 



> 



2g(7-^) 
7 + g 



> 1 



(5) 
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Fig. 4. (a) Phase diagram showing the isolated islands in the model's parameter 
space where it is possible to tunnel from mortality into immortality in the case 
a = 0. The vertical lines denote 7 values for: UK, Sweden, Australia (A); US, 
France, Italy, Canada, India (B); Spain, Greece, Romania (C). For E/C = 1.5 all 
the above countries lie on the first island, (b) Time evolution of wealth distributions 
for the case of life extension with tunneling (7 = 2, a = 0, C = 2, = 3). In this case 
the introduction of life extension leads to the emergence of a bi-modal distribution 
of wealth. 

These inequalities specify a series of isolated islands within the (7, E/C) 
plane in which it is possible to tunnel from mortality into immortality. These 
islands are indexed by integer n. For the simplest case a = they are the 
areas between curves ^ = ("+1)7 g^^^ e _ — n'~i{'y-i) ^ r^j^^ critical 

G 2(7— a)— n— 1 C 2(7— q)(7— 1)— 717 

points 7„ = (1 + \Jl + 2n{n + l))/2 denote the start of the n-th island. It 
follows that mortals and immortals necessarily form segregated populations if 
7 is smaller than the golden mean 71 = (1 + \/b)/2. 

Fig. shows the first few islands in which tunneling is allowed for the case 
a = (similar graphs follow for other values of a). The vertical lines in Fig. Hti 
denote 7 values for: UK, Sweden, Australia (A); US, France, Italy, Canada, 
India (B); Spain, Greece, Romania (C). These have been calculated using data 
for adjusted annual growth of these countries from 1960 to 2000 and assuming 
that one generation in our model corresponds to 30 years p^. For E /C = 1.5 
all the above countries lie on the first island in the above phase diagram. 
The much higher growths of China and Singapore are also consistent with 
tunneling between mortals and immortals for the same value of E/C since 
they lie on the second island. Fig. |Dd shows that life extension profoundly 
influences the distribution of wealth. For the society with 7 = 2, a = 0, 
C = 2, and E = 3, tunneling takes an initial Gaussian wealth distribution 
into a bi-modal one. Wealth distributions of this type are very similar to those 
of existing highly segregated societies in which life expectancies at birth differ 
significantly between the rich and the poor [31j. This is an indication that 
the presented model, although substantially simplified, captures key aspects 
of realistic processes. 
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Fig. 5. Profit of pharmaceutical companies as a function of the unit price E. The 
society shown has 7 = 2,a = 0,C = 2. The initial wealth distribution was a 
Gaussian centered at m = 2 with width a = 1. The t = 1, 2,4 time slices are shown 
from bottom to top. The maximal profit for the producers of life extension is for 
£; PS 3, i.e. for E/C ss 1.5. 

Similar changes may be found when looking at fertility rates of mortals and 
immortals, as well as the overall size of the population. In contradistinction 
to what one might naively expect, the introduction of life extension does not 
speed up population growth. In fact, for realistic values of 7 the size of the 
population generally stabilizes. This gives us interesting examples of societies 
with sustained economic growth but without a spiraling population explosion. 
It is not difficult to see that this is a consequence of the dynamical equilibrium 
between two phases in the model. In fact, this uncovered non-trivial behavior 
within a simplified model is the essence of how physics can contribute to our 
understanding of economy and society in general. Namely, effective models in 
physics (e.g. the Ising model) are of value not because they encode the detailed 
phenomenology, but because they capture key qualitative relations between 
dynamical quantities like the one above, providing insight needed for deeper 
understanding of the underlying phenomenology. 

So far we have looked at life extension from the consumer's perspective. We 
now briefiy look at the profits of the pharmaceutical companies selling the 
life extension product. Each individual purchase of life extension increases 
the profit of the life extension companies by E. We assume here that all 
the R&D expenses of developing the product have already been covered and 
that the actual cost of manufacturing the product is negligible. Summing the 
purchases over the whole society we get the time dependence of the total profit. 
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An example of this is illustrated in Fig. [5] for a society with 7 = 2, a = 0, 
and C = 2. From the figure we see that the maximal profit determines that 
E/C ~ 1.5. Note that this is also the horizontal line in the phase diagram 
in Fig. S^. We see, therefore, that the economic requirement of maximizing 
profit of pharmaceutical companies is consistent with the political requirement 
of easing social tension through de-segregating mortals and immortals, i.e. 
through allowing tunneling into immortality. 



4 Concluding remarks 

We have presented and solved a simplified model that analyzes the con- 
sequences of (repeatable) life extension on fertility, population growth and 
wealth distribution. When life extension is absent the model correctly repro- 
duces observed time dependence of wealth distributions, and abrupt declines 
in fertility rates. We have analyzed in detail the introduction of life extension 
to the model and have found it to be a novel commodity which profoundly 
influences key aspects of society. Of particular interest is the emergence of two 
distinct phases: societies in which mortals and immortals are segregated, and 
societies in which economic factors allow descendants of mortals to "tunnel" 
into immortality. 

The analysis of simplified models such as ours is but a first step in a pro- 
cess that could ultimately help in forming important future policies, e.g. those 
to do with the pricing of pharmaceutical and medical products and services, 
wider health-care and insurance policies, etc. As is well known, these issues can 
have profound effects on the stability of societies and their economic growth, 
and have for this reason attracted much attention. An important recent ex- 
ample is the decision of the Brazilian government to bypass the copyright on 
US-produced AIDS drugs [34] in order to be able to treat significant part 
of the AIDS-infected population, and to avoid political instability that may 
arise from this problem. Models such as ours have the possibility of leading 
to rationally thought-out policies, allowing society to make critical choices ac- 
ceptable to its members. However, in order to do this they necessarily need 
to be followed up by the development and analysis of a series of richer mod- 
els incorporating more realistic assumptions. We have already commented on 
some directions in which this process of model building needs to go when 
we discussed the assumptions within our basic model without life extension. 
The introduction of life extension will further affect matters such as depen- 
dency ratios, or the effects of overlapping generations. The issue of work-leisure 
tradeoffs can also play an important role in the dynamics of a society with 
life extension. On the other hand, the very introduction of life extension could 
greatly influence our attitudes towards work and leisure in prolonged life spans. 
Also, the decisions based on risk assessment will necessarily undergo a qual- 
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itative change when made from longer time perspective offered by repeated 
hfe extension. Social and economic strategies will also change accordingly. 

In the presented model the economic growth, parental investment, and price 
of life extension were all externally determined. More realistic models should 
include the boot-strap influence of population growth and wealth distribution 
on these parameters. We plan to pursue this generalization in a future publi- 
cation. Another interesting extension of the model would be to consider the 
interaction and co-existence of two parts of society having different parental 
investments. 
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